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Початкове наближення для обчислення складного щільнісного контакту 


Запропоновані аналітичні моделі гравітаційного поля та горизонтально-шаруватого геологічного сере- 
довища з кількома густинними межами розділу шарів під поверхнею. Модель поля отримана з рівняння сили 
тяжіння, що випливає з розгляду плоского нормального потенціалу сили тяжіння у локальному околі точки 
вимірювань. Введено нову модель середовища - клас Чорного контактних поверхонь. Подано нову схему ви- 
значення початкового наближення густинного контакту і його асимптот у багатошаровому геологічному 
середовищі за допомогою кількох 'контактних'алгоритмів. 
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Основний напрям наукового прогресу - створення 
загальних математичних моделей в геології, геофізи- 
ці і надалі - створення геолого-математических тео- 
рій... (Елланський М.М., 2003, ммим.ттеїї.пагод.ги) 


Вступ. Нині чималі зусилля геофізиків спрямовані на переінтерпретацію картографічного матеріалу, от- 
риманого за даними зйомок 30-річної давності, зважаючи на різке обмеження обсягу геологорозвідки. Про- 
те, через підвищені вимоги до точності і ефективності інтерпретації, така переобробка великих обсягів гео- 
фізичних даних потребує нових підходів до чисельного обгрунтування методів, використаних при підготов- 
ці і обробці подібного картографічного матеріалу. Цілком логічно об'єднати розвинені на попередньому 
етапі розвитку геофізики математичні конструкції із сучасним інтерактивним програмним забезпеченням, 
щоб комплексно вирішити проблеми оцифровування, візуалізації і тлумачення геофізичних даних. Пробле- 
му тлумачення даних, найважливіша для застосувань (втім, тісно переплетена на рівні первинної обробки 
із першими двома), вирішують останнім часом у руслі перебудови аналітичної і обчислювальної бази теорії 
інтерпретації потенціальних полів, яка забезпечує становлення нової парадигми в теорії і практиці інтер- 
претації. 

Сучасні методи обробки геофізичних даних розвиваються із прицілом на створення цифрових аналі- 
тичних моделей поля і геологічного середовища і високоточних алгоритмів, орієнтованих на застосуван- 
ня в пакетах ГІС, особливо при відображенні вхідної інформації і результатів її трансформації та зістав- 
ленні із апріорною інформацією, також поданою в графічному виді (т. зв. візуальна інтерпретація, рис. 1). 
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Рис. 1. Карта гравітаційних аномалій України в середовищі вільно поширюваної ГИС ОСиапіит бІ8. 


Власне, "пошарова" обробка являє новий етап якісної інтерпретації, в якому можливості подання ін- 
формації у вигляді електронного журналу вимірів (13) можна використати для попередньої обробки ін- 
формації в середовищі спеціалізованих пакетів загального і спеціального призначення. Як приклад такої 
попередньої обробки, можна вказати процедури: виділення регіонального фону, наприклад, шляхом міні- 
мізації нев'язки між набором даних і деяким поліномом в середовищі Маніаб; попередня фільтрація вимі- 
рів в середовищі Ехсеї; знаходження особливих точок поля методом зворотної інтерполяції полінома Ла- 
гранжа в Маріе, редукції поля за допомогою алгоритмів В.М. Страхова, реалізовані в Руїпоп, і т.п. Одним 
із найбільш вдалих прикладів цифрової якісної інтерпретації можна вважати приклад обробки значень 
поля сили тяжіння в зв'язці ПЗ Маріпіо - СогеїОгам (рис. 2) | 131, що в жодному разі не виключає їх заміни 
в осяжному майбутньому на ореп-8оигсе аналоги. 
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М просмотр резутьтетов 
Рис. 2. Попередня обробка гравітаційних аномалій в зв'язці Маріпіо - СогеїОгам. 

Проте ми торкнемося іншого напряму у рамках детермінованого підходу до інтерпретації даних поте- 
нціальних полів - розробки аналітичних конструкцій геофізичного поля і геологічного середовища. Щодо 
математичних конструкцій поля, відсилаємо зацікавлених читачів до праці (12), а серед моделей середо- 
вища автору близькі моделі густини типу контакту, одну з яких пропонуємо вашій увазі. 

В розробці таких конструкцій заради підвищення точності інтерпретації даних геофізики обов'язково 
слід врахувати умови адекватності постановок сучасним вимогам до отримання вихідних геофізичних 
даних і їх обробки. В цілому, суть цих вимог така: 

. висока точність і технологічність алгоритмів і програм, їх переносимість; 

. адаптація моделей поля до задання початкових вимірів на коротких профілях; 

. адаптація моделей середовища до його складного (нелінійного неоднорідного) стану. 

Визначення контакту у відомих моделях середовища Страхова, Нумерова, Маловичка (| 1) замкнена 
на обчислення різниці ундуляцій деякої теоретичної моделі середовища від деякого опорного її елемен- 
та, від якого ці ундуляції не надто сильно ухиляються. Каверза полягає в тому, що при цьому необхідно 
заздалегідь знати деякі елементи геометрії (глибину, форму, орієнтацію) цього опорного елемента (мо- 
дельного тіла або контакту), як початкового наближення ітерацій. Від вибору вдалого початкового набли- 
ження залежить, в цілому, ефективність будь-якого із існуючих ітераційних алгоритмів (71. 

Ми вивели для горизонтально-шаруватого середовища з кількома густинними межами шарів аналіти- 
чні моделі поля і геологічного середовища, в яких ці опорні елементи розраховуються в процесі ітерацій. 
Модель поля є рівнянням сили тяжіння, і отримане з розгляду нормального потенціалу сили тяжіння в 
локальному околі точки вимірів в двовимірному поданні. Розміри локального околу істотно залежать від 
розмірів тяжіючої області і точності апроксимації нормальної складової потенціалу |51. 


Аналітична модель поля. Нехай у тривимірному необмеженому Евклідовому просторі ЕФ) дві області 
Сі С, вкладені одна в одну і відстань між їх межами дС; і дС» - 4ї5(ОС,,дС)20 при б, б, «б». В 
області С» початок координат і точки області 8" «6,26», С» - б» одС», позначмо через С - (а |, 
а точки замкнутої необмеженої області С'«ЕФЛС, дб'з дС - через ха (ххх), тоді 

5 у2 
або) фе (Ми о, збо 
ізі 
Нехай область С» заповнена масами з щільністю дб» (2), а область С - С, ХС» - масами з щільністю 


5.(2) (рис. 3). Позначимо через п (х) одиничну нормаль до поверхні ОС" , зовнішню щодо С; тоді нор- 
мальна складова сили тяжіння дорівнює 
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де 5(2) - 65(2)- 52), Є є б», (2)- 0, Є є Су, /- гравітаційна константа. 
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Рис. 3. Модель контактного середовища в класі Чорного. 
Взагалі кажучи, силу тяжіння описує функція 


цОх) - |втад у (х) - / Ї (є )втад ооо е ДИ, (2) 
р 
оскільки значення сили тяжіння - це значення модуля градієнта потенціалу |41, а не складові тяжіння, у то- 
му числі і (1). 
Спростимо вираз (1), розглянувши малий локальний окіл 0 точки х. У цьому околі, з відомою долею 


наближення | |51, ділянку СС. СОсдС" межі ОС" вважаємо "необмеженою" площиною: 
ОС з їх хат РА РОВЕН з соп5і, (хіьх») є З де 5, - 0 0.дС". Розташувавши в цій площині прямокутну 
систему координат так, щоб вісь охз була перпендикулярна до площини, отримаємо 

я) є є, хе (о, хазха), бо (б,боьба), Оз ОпеТ, 0 заое, дс Од. 
Нехай в околі 07 функції (2) сопяї, Є є С, с 07, 65(2)- сопіі, СЕ О оС, б-бь -дбі,а С» - влас- 
на підмножина 0 . Фактичні розміри області 0 залежать, передусім, від розмірів тяжіючої області С, , що 
вивчається, і точності апроксимації виразу |51, тобто 
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Якщо є»0 досить мале, то СО вибрана вдало. Межу С вибираємо досить гладкою, для якої 
дС» з Їх же хор) С рхи)є 05, 1,7,К - 1,2,3; 05; « дС» ух; є сб) (д5, ). 
За таких припущень про середовище отримуємо потенціал тяжіння 





цаоз)є піз -к пооздс , Ка є (4) 
3 ог ОХ3 ох)? 
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Якщо маси, зосереджені в СО , генерують аномалії и(х,)), витягнуті, скажімо, уздовж осі ох» , їх краще 
вимірювати по профілях, паралельних од; . Формула (4) малопридатна для обробки таких вимірів. Проїін- 
тегруємо вираз (4) по б» і після ряду аналітичних перетворень отримаємо: 


сх 93 -Хз 
м) "г (2, 7-х 5 "В (а 7 Хз у зо 


де 5 - проекція б, на площину охухз, К.- Іс. При ху 20, ба 52), ас «5 отримаємо модель 
контакту у вигляді 
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де 209 з З) і 20) п ЕХ), а «2 «Р - криві, що обмежують знизу і згори тяжіючу плоску область 5. 








Аналітична модель середовища. Нехай С  - горизонтальна смуга, обмежена згори і знизу прямими 
де, сей, 04 «А «оо, яка розділена на 2 шари із щільністю д; і б, нескінченно гладкою кри- 


вою ОС", зірчастою щодо нескінченно віддаленої точки, з асимптотою 2-й, й, «п « й,, причому 
«(29 -| є а»д, |», «)єдс (6) 


Якщо значення Л і величина стрибка щільності Дб з б» - б; відомі, то межа СС розділу шарів за зовніш- 


нім полем відновлюється однозначно. Розв'язок оберненої задачі про контактну поверхню з класу (б) єди- 
ний |81, якщо відомі стрибок щільності Дд і скінченна дуга кривої 2(2) з (6). 





Введемо спеціальне позначення для класу (6) о ЛЬ (-20, со) і назвемо його класом Чорного - на честь 


вітчизняного науковця-геофізика А.В. Чорного, який заклав підвалини чіткого алгоритму дослідження 
функціональних множин в царині контактних поверхонь. Очевидно, що 


СНО РУ (оо, со) є ЗР (со, оо) , 8 а, Б) с 5190, оо). (7) 
У такій постановці важливою є проблема розділення полів, яка є самостійною задачею і виходить за ра- 
мки цієї статті. 
Щоб розв'язати нелінійне інтегральне рівняння (5), яке описує контакт, обмежений асимптотами у плос- 
кій області С» , його слід тлумачити як функціонал щодо шуканої функції хз(21) на деякому лінеалі. Для 
цього функціоналу нескладно знайти приріст Де(х, 2) - я(х, 2-- й) - є(Х, 2), отримуючи 
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Відновлення контакту. На основі (8) можна збудувати кілька способів відновлення контакту хз(сі), 





ассє«Ь,зазаданим полем є(х), сххх«а, |а, Р| с |с, 41, пале5(а, Б| « плез|с, а). 
Спосіб Ньютона. Вважаючи відомими а,?, Н , в результаті розкладу (14) в ряд по хз(С|), та подаль- 
ших перетворень знаходимо АД (2) є 2,ц(4)- 7, а 
Б 22 (2) -ху-нН? 21 и ) 
Ре 5 ає-ксоз| о : 
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Альтернатива. Якщо вираз (8) записати точністю до і? (2) та яз (2) , то після ряду перетворень от- 
римаємо, відповідно, способи "контакт 1"та "контакт 2" 
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Рівняння (9-11) можна вирішувати ітеративно, якщо задано 25(2): 21(5) З 20(4) 4 Йо(5). 


Проте можна одночасно визначати не лише контакт 2 - (2), ас хь, але і"асимптоту" Н - соп5і, 
модифікувавши дещо вказані вище способи. 
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ав -х я зно о на ФІЇ 
Початкове наближення контакту. Нехай крива є(х), ссхоха маєна |с, | кілька екстремумів у то- 
чках с, є|с, д|, і - 1,2,..., п. Оскільки ми поклали б » 0, то кожному екстремуму відповідає однойменний 


екстремум 2-2(2), асо ар. Визначимо в множині (ехг я(х)) найбільший і найменший екстремуми 
Іс, 4) 


(найбільший максимум 2; і найменший мінімум 2; ), тобто 
85 - 8(с;) - 8аріемг 800), 8; з 8(сз) - іп ехіг 8 (х)). 
Гс, 4) Гс, 4) 
Позначимо через 25 з 2(с;) - 5ар(ехг с(2)), 2; з 2(с,) з ШЕ ехіг с(2)) відповідні екстремуми с. Тут абс- 
Іс, а) с, 4) 


циси 88 і 25 (8; і 2; ) співпадають, що не завжди вірно, але в даному разі так. Початкове наближення кон- 
такту можна задавати, щонайменше, двома способами (02 спрямована вгору) 


25 2 20(4)2-Но 5 дє(5, Но), ассхр, (15) 
29 З 20(8) 7 -Но З Коде(с), акохь, (16) 
де Но»0 - константа, що задається, - "асимптота" функції 2, на лівому кінці проміжку (а, Р), 


ДЕ(О з 2(2)- 8(а), Де(2, Но) - аналітично продовжене значення Де(2) на рівні -Но, Ко - К(Но) - 
коефіцієнт, що підлягає визначенню. 
Визначення Ко у підсумку зводиться до визначення початкових 2; і с; призаданих Н, 85 Ї 8;: 
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з:А(-с ни Р тб; іно, - 2Н мав оту онов (18) 
Ш-ер'яар |" сіна зоре) зі - с; (1-с;) 
Утворимо функцію /(2) -(1- с) ІА - су? ча Інесса) 22 агоів рдеруляме 
-с Га 
а-стіа-с?нНн?)-сінес -Н2)-2Н нти ро. (19) 
НІ -с(1- с) 
Наближення 2, якогось із значень 25 і с; визначаємо за схемою Чебишева (11: 
" 2 
аазо Год Годо, (20) 





Га) | 27") 
Процес (20) при вдалому виборі 2, швидко збігається: 3-4 ітерації для отримання 2 
що узгоджується з результатами, отриманими в (1). 


0) з точністю до 107, 


Обчисливши таким чином для кожного заданого значення Н;, і - 1,2,..., п, лівої "асимптоти" шуканої 


кривої 2 - (2), ас 2 сЬ, значення параметрів 29 , а, НІ, і- 1,2,..., п, визначимо 


їз 
и? з10) - НІ, цу є20)-НІ, АН, еНІ -НІ 

58 - 8(Р)- 8(а), бд8пих 7 85 3 8(), б8 іп й 

| (і) 

и 








сіно Зк ау ЇХ (і) (6) (і) що 
і знайдемо шукані коефіцієнти Ку - -УХК; ,деки- Мо 


е їі |б в пах 








Відшукання початкових наближень на конкретних моделях контактів буде матеріалом наступної статті. 


Висновки. Введено в науковий обіг нові аналітичні конструкції для обчислення складних контактів у ви- 
гляді виразів (17), (18). Для їх обчислення запропоновано процес простих ітерацій (20). При чисельному 
моделюванні вказані вище конструкції виявляють кращу збіжність при обчисленні за способом Ньютона 
(20) порівняно із класичними прикладами |З), але їх практична ефективність вимагає перевірки на підго- 
товленому матеріалі польових гравіметрових знімань.. 
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